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ABSTRACT 

Aims. DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition 

and energetics of the circumstellar material during this pivotal stage. 

Methods. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 /im as part 

of the DIGIT Key Program. 

Results. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO 

from 7=14-13 to 38-37 (EJk = 4080 K), water from levels as excited as Jk^^k^^ = 7o7 (EJk = 843 K), and OH lines up to EJk = 290 K were 

detected. 

Conclusions. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk 

of 0.03 Mq, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, 

atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but 

some emission could be due to C-shocks in the walls of the outflow cavity. 



1. Introduction 

One of the pivotal points in the formation of a star occurs when 
the star and disk system emerges from the collapsing envelope, 
signaling the end of the main accretion of stellar mass. The tran- 
sition from Stage I protostar, embedded in a substantial enve- 
lope, to a Stage II configuration, with only a star and disly is 
poorly understood. 

IRAS 12496-7650 is the brightest infrared (IR) and sub- 
millimeter (submm) source in the Chamaeleo n II star-formin g 
region, at a distance of 178 + 18 pc ( Whi ttet et aH Il997h . 
Follow-up studies to the discovery of IRAS 12496-7650 by 



iHughes et al.' (11989^ 1991 ), along wit h subsequent observations 
(e.g..lWh ittet et al.l ll997llYounget al . 2005, P orras etal.1 120071 
ISpezzi et al.l l2008). identified the variable Herbig Ae star DK 
Cha as the optical counterpart. The central star mass is esti- 
mated at 2-3 Mo (ISpezzi e t al.|2008D, with a bolometric luminos- 
ity betw een 24 (Sp ezzi et al.l2008l) and 35 Lq ([van Kemp en et al.l 
l2009ah . The star is surrounded by a molecu lar outflow and pro- 
tostellar envelope (Ivan Kempen et al] 120061: l2009ai) . The spec- 
tral energy distribution (SED) is roughly flat in wavelength 
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* Herschel is an ES A space observatory with science instruments pro- 
vided by the European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

' See, e.g.. iRobitaille et al] ( l2007h for a definition of the Stages 



from near-infrared to mid-infrared dHughes et al.lll989l) . Using 
the Tboi definition and including longer wavele ngth data, (e.g. 
iHenning et all I1993L Ivan Kempen et alJ l2009ah . DK Cha has 
Tboi = 580 K, making it a Class I source, but near the Class 
II boundary at 650K. Due to a n approximately face-on v iew- 
ing angle down the outflow cone (Ivan Kempen et al]l2009ah : one 
can observe the protostellar disk directly, unobstructed by the 
surrounding envelope. This viewing angle also allows the op- 
tical emission to escape. The observational SED Classes are 
often assumed to coiTespond to the physical Stages, but this 
is by no means obvious when aspherical mo dels are exam- 
ined (iRobitaineitaD l2007l ICrapsi et al.1 l2008l IPu nham et aL 
2010") or tracers of dense gas are considered (van Kemp en et al. 
2009b). The same physical object can have different SED clas- 
sifications depending on viewing angle. Thus, it is interesting to 
study an object that is near an SED classification boundary and 
that also appears to be transitional in physical stage. 

Stron g far-IR continuum a nd hne emission have been 
detected (Giannini et al.' '1999', van Kempen et al. '2006'). As 
discussed by van Dishoeck (2004), Nisini et al. (2005), and 
Ivan Kempen et al.l (12(^1 0) for the case of the similar source 
HH46, far-IR atomic, ionic, and molecular lines can constrain 
the properties of the disk, envelope, and outflow. CO emission 
lines between 7=14-13 and 7=19-18, as well as [Oi] and [Cii], 
were observed from DK Cha with ISO-LWS, but other possi- 
ble molecular species, such as OH and H2 O, were undetected 
dLorenzetti et alJll999llG"iannini et al.lll999l) . 
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Table 1. Detected gas lines within the PACS spectrum of DK 
Cha. 



Species 


Transition 


A" 


E„ 


FWHM 


Flux* 






fim 


K 


jjm 


Wm-2 


[Oi] 


'Pl-'P2 


63.18 


227.72 


0.08 


3657.7 


OH 3/2-3/2 


9/2- - 7/2+ 


65.13 


512.1 


0.05 


56.0^ 


OH 3/2-3/2 


9/2+ - 7/2- 


65.28 


510.9 


0.07 


52.32 


CO 


38-37 


69.07 


4079.98 


0.06 


42. l2 


CO 


37-36 


70.91 


3871.88 


0.09 


59.6^ 


o-HjO' 


7o7-6i6 


71.95 


843.5 


0.06 


37.12 


CO 


36-35 


72.84 


3668.96 


0.06 


40. 1^ 


CO 


35-34 


74.89 


3471.44 


0.1 


57.0^ 


0-H2O 


423-312 


78.74 


432.2 


0.04 


19.8^ 


OH 1/2-3/2 


1/2- - 3/2+ 


79.12 


181.7 


0.06 


282 


CO 


33-32 


79.36 


3092.45 


0.03 


17.52 


CO 


32-31 


81.81 


2911.15 


0.05 


126.2 


0-H20 


6i6-5o5 


82.09 


643.5 


0.03 


96.2 


CO 


31-30 


84.41 


2735.28 


0.09 


336.0 


OH 3/2-3/2 


7/2+ - 5/2- 


84.42' 


291.2 


- 


- 


OH 3/2-3/2 


7/2- - 5/2+ 


84.60 


290.5 


0.09 


301.4 


CO 


30-29 


87.19 


2564.94 


0.09 


317.5 


CO 


29-28 


90.16 


2399.93 


0.14 


545.8 


CO 


28-27 


93.35 


2240.30 


0.07 


261.4 


CO 


27-26 


96.77 


2086.22 


0.10 


330.1 


CO 


25-24 


104.44 


1794.23 


0.2 


314.5 


0-H20 


221-I10 


108.08 


194.1 


0.15 


126.9 


CO 


24-23 


108.76 


1656.55 


0.15 


260.0 


CO 


23-22 


113.46 


1524.19 


0.2 


477.4 


0-H20 


4l4-3o3 


113.54' 


323.5 


- 


- 


CO 


22-21 


118.58 


1397.43 


0.17 


285.2 


OH 3/2-3/2 


5/2- - 3/2+ 


119.23 


120.7 


0.18 


<40 


OH 3/2-3/2 


5/2+ - 3/2- 


119.44 


120.5 


0.17 


<72.3 


[Nil] 


'P2-3P, 


121.90 


188.20 


0.12 


45.7 


CO 


21-20 


124.19 


1276.10 


0.17 


357.1 


P-H20 


404-3l3 


125.35 


319.5 


0.12 


42.8 


CO 


20-19 


130.37 


1160.24 


0.16 


319.0 


CO 


19-18 


137.20 


1049.89 


0.18 


338.7 


P-H20 


3l3-202 


138.53 


204.7 


0.17 


103.0 


CO 


18-17 


144.78 


945.01 


0.15 


350.4 


[Oi] 


^Po-'P2 


145.53 


326.58 


0.16 


146.7 


CO 


17-16 


153.27 


845.63 


0.15 


386.0 


CO 


16-15 


162.81 


751.76 


0.18 


466.5 


OH 1/2-1/2 


3/2+ - 1/2- 


163.12 


270.2 


0.24 


91.0 


OH 1/2-1/2 


3/2+ - 1/2- 


163.40 


269.8 


0.19 


90.0 


CO 


15-14 


173.63 


663.37 


0.16 


413.0 


0-H20 


303-2l2 


174.63 


196.8 


0.19 


208.0 


0-H20 


212-101 


179.53 


114.4 


0.11 


31.0 


0-H20 


221-212 


180.49 


194.1 


0.12 


34.9 


CO 


14-13 


186.00 


580.51 


0.20 


437.2 



" Laboratory wavelength. * Flux summed over all pixels in 10""* W m"2 Typical 
line flux uncertainties are dominated by calibration uncertainties (~ 30-50%) and 
fringing (0.1 Jy at the longest wavelengths to 0.6 Jy at the shortest). 
' Lines blended with above line. - Line detected only in central spaxel and should 
not be compai'ed to the fluxes of the summed spaxels. ' Doublet unresolved. 



The large beam of ISO-LWS (90") and its low spectral re- 
solving power (A/AA = 200) were insufficient to determine the 
origin (e.g., envelope versus outflow) of the CO emission. The 
preferred picture was of a small, quiescent hot core of ~300 AU 
around the protostar (iGiannini et al.lll999h . but recent observa- 
tions of CO 4-3 and 7-6 indicate that this warm CO emission 
either originates outside the inner few thousand AU or cannot be 
heat ed simply by r adiation from the central star (passive heat- 
ing) (I van Kempen et al...2006.) : the far-IR and submm lines have 
different energies and critical densities, so they may originate in 
different regions. 

The Key Programme, Dust, Ice, and Gas in Time (DIGIT) 
aims to study the evolution of these three constituents through 
the star formation process by observing objects in different 
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Fig. 1. SED of DK Cha, using WFI I-band data, 2MASS (J,H,K), 
Spitzer-IRAC (3.6 and 4.5 fim. The other bands, 5.8 and 8.0 
jt/m are not shown due to saturation), Spitzer-MIPS (24, 70 
and 160 fim), APEX-Laboca (870 fim) and SEST (1.2 a nd 1.3 
mm). Overplotted are the best-fit model from the Robitaill e et al.l 
(12007 ) grid and the PACS and Spitzer-IRS spectra. The MIPS 
70 fim is corrected for saturation; the 160 fim for large-scale 
emission. The solid line is the full model, while the dashed line 
shows the envelope contribution, the dot-dashed the disk con- 
tribution, and the dotted the star. Photometric points are shown 
as diamonds. The short-wavelength points are not well matched. 
Fluxes at short wavelength (<3 fim) are variable. 



stages. This pape r, along with a com panion paper on HD 100546, 
a Stage II object (ISturm et al .'^2010). present the first results from 
the DIGIT Science Demonstration observations. 



2. Observations 



The launch of the Herschel Space Observatory (|PilbratJl2010l) . 
with the Photodetector Array Camera and Spectrometer (PACS) 
instrument (Poglitsch 2010), enables far-IR spectra with high 
spectral resolving power and good spatial (~9") resolution. The 
range scan/SED mode provides a full spectrum (55-210 yum). 

DK Cha was observed for a total of ~4 hours on 10 
December 2009 by the PACS spectrometer on the Herschel 
Space Observatory in pointed range-scan spectroscopy mode 
(obsid 1342188039 and 1342188040). An image slicer rear- 
ranges the light into a 5 x 5 array, with each spatial pixel (spaxel) 
covering 9'.'4x9'.'4. A grating disperses the light onto two 16 x 25 
pixel Ge:Ga detectors. Observations were taken using nodding, 
in which two off-positions separated by ±6 arcmin in elevation 
were taken as 'blank' sky, from which the signal is subtracted. 
Consequently, spatially diffuse lines (e.g., [Cii]) are subtracted 
out. 

The effective spectral resolving power was about lOOO-i-/- 
500, lower than was expected. This is due to the pointing er- 
ror, which is perpendicular to the slit, as well as due to fringing 
issues on the continuum. The point-spread function of PACS in- 
creases towards longer wavelengths and is larger than a single 
spaxel, but is corrected in post-pr ocessing using the PACS spec- 
trum of the pointlike HD 100546 (iSturm et alJlJOlOh to calculate 
the extraction efficiencies. 

The data were reduced using the Herschel Interactive 
Processing Environment (HIPE v. 2.4). The observation was 
slightly mispointed (~5-6") making it hard to assess whether 
emission is extended and causing wavelength shifts in some 
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Fig. 2. The spectrum from 55-190 fim of all spaxels added together, after removal of the continuum and rebinning to 0.5 fim 
resolution. Several dips or bumps (e.g., around 75 jum) are instrument artifacts. Inserts show expanded views of parts of the spectrum, 
with line identifications. 



spectral ranges. Because of the mispointing, uncertainty in the 
Spectral Response Function, light leakage problems, and uncer- 
tainties in the wavelength dependence of the PSF, the spectrum 
was trimmed shortward of 56 yum, longwardof 190jum, and from 
95-101 fim. Two spectra were used to identify lines. One is from 
the central spaxel alone. For a second spectrum, used for most 
lines, fluxes from all spatial pixels were summed to correct for 
the PSF and mispointing. We will focus on gas lines 



Table 2. Best-fit from Robitaille models 



DK Cha characteristics 



Envelope Mass 


0.03 Mo 


Al^ 


7.9 


Disk Mass 


0.03 Mo 


Inclination 


< 18° 


Luminosity 


29.4 Lo 


Disk Inner Radius 


0.3 AU 



' Total Ay of which Av=5.2 mag. is attributed to the envelope at an 
inclination of 18° and 2.7 mag. to more extended cloud. 



3. Results 

3.1. Continuum 

The total Spectral Energy Distri bution (SEP) in Fig. 1 includes 
photome try from the lite rature (lAlcala et al.ll2008l ISpezzi et aTl 



[2008, E vans et al.l l2009l Nefs et al., in prep), combined with 
Spitzer-IRS spectroscopy and the PACS results. The photometric 
points were fitted by the SED fitter routine from .Robitaille et al.. 
(l2007h . with results for parameters in Table|2] Due to the extinc- 
tion and variability at short wavelengths, only points longwards 
of 3 fim were taken as input for the fit.The best-fit model is over- 
plotted and split into the contributions from the star, disk, and 
envelope. The short- wavelength IRS continuum is dominated by 
the flux from the disk, the longer wavelength IRS continuum 
covers the transition from disk to envelope, and the PACS con- 
tinuum arises in the envelope. The model underestimates the ob- 
servations in the near-lR and mid-lR, indicating that the disk 
emission is underestimated in the model. 

The best-fit models of the SED grid indeed confirm the as- 
sumption of a nearly face-on configuration <18°, with relatively 
low and equal envelope and disk masses of 0.0 3 Mp, indicative 
of a source in transition from Stage 1 to Stage II dRobitaille et al.l 
12007 ). The PACS continuum was not used to constrain the fit; 
the consistency with the predicted model (differences of <10%) 
is thus encouraging, especially considering the calibration un- 
certainty of 30-50%. 



3.2. Line Inventory 

Table[T]lists all gas lines that have been identified. Fig. 2 shows 
the spectrum after subtraction of a polynomial fit to the contin- 
uum. We us ed a combina tion of the ISO line lisQ and the list 
published in iLerate et al.l ([^06) for the Orion bar More lines 
may be detectable in this spectrum after PSF deconvolution. 
Some lines in Table[T]are only detected in the central spaxel. Due 
to the fringing in other spaxels at short wavelengths and the PSF 
problems discussed above, these do not show up in the summed 
result of Fig. 2. Using spectra both from the summation of all 25 
spaxels and from only the central spaxel, a large range of CO, 
H2O and OH lines were detected, ranging in excitation energies 
from a few hundred kelvin to over 3000 K. In addition, [O i] and 
[Nil] were identified. Emissio n from [Cn] at 157.7 /im is not 
seen. That line was detected bv lLorenzetti et al] ( 1199 9), but was 
believed to be extended emission, associated with the galactic 
ISM. This component is subtracted out in our observations due 
to the nodding, which was not used in ISO-LWS. 

3.3. Gas analysis 

Far-infrared CO lines were previously detected bv lGiannini et al.l 
(1999), but only up to CO 7=19-18. Assuming flie Unes to be 
optically thin, those could be fitted with temperatures ranging 
from r = 200 K to r = 750 K and densities from n = 4 x 10^ 



http://www.mpe-garching.mpg.de/iso/linelists/Molecular.html 
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Fig. 3. Upper: CO line fluxes plotted versus /up- The ISO line 
fluxes are shown with triangles. A clear excess can be seen in the 
CO line fluxes from /up -2 onwards compared t o the warmest 
model of 750 K proposed bv lGiannini et alJd 19991) (dashed line). 
Lower: CO Rotational diagram, showing two distinct popula- 
tions with a break around 1500 K in energy. The blended lines 
are excluded from the fit and the plot. CO lines with 7up larger 
than 3 1 are excluded from both plots, as these lines are detected 
only in the central spaxel. 



5 X 10^ cm"^. The high T, low n model, which 
produces the most emission in the high-7 lines, is overplotted in 
Fig. 3. As can be seen, this model severely underproduces the 
CO lines from higher transitions seen by Herschel. The other 
proposed models produce even less high-/ emission. 

The CO level populations up to 7 = 31 (Fig. 3), can be fitted 
with two distinct rotational temperatures {Troi), with the break at 
about iiup = 1500 K. Most of the mid-J CO lines can be fitted by 
Troi ~ 380 K, while high er excitation lines are be st fitted with 
Tr„, ~ 1 365 K. In HH46 (Ivan Kempen et al.l l2010l). UV photon 
heating (ISpaans et al.ll 19951) is dominant for lines with energies 
below 1700 K, while a non-dissociative C-shock dominates the 
emission from higher excitation lines. A similar model is plausi- 
ble for DK Cha, and the issue will be explored in a forthcoming 
paper (Green et al. in prep). 

Water was detected in a wide variety of transitions, up to 
7o7-6i6 at 71. 93 ^um with (E/k - 8 43 K). Water was not detected 
in ISO-LWS (JGiannini et al. ' 1999) because the lines are 5 times 
weaker than the CO lines. The discovery of the 7o7-6i6 fine in- 
dicates high temperature, as confirmed by the high-7 CO Unes. 
Given the possible identification of shocked CO and the geome- 
try of DK Cha, it is likely that the water emission arises at least 
in part from either a non-dissociative C-shock or UV-heating of 
the quiescent material along the outflow cavity walls. 

In addition, at least 9 lines of OH are detected. Since oxygen 
is rapidly transformed into H2O at temperatures higher than 100 
K and frozen out below that, either a very strong UV field or a 
dissociative shock is needed. OH is a m ajor coolant for dissocia- 
tive shocks dNeufeld & Dalgarnoll989i) . This dissociative shock 



is accompanied by a C-shock, which may explain the CO and 
H2O analysis. Models from Wampfler et al. (in prep) also sug- 
gest that C-shock densities are often too low to both account for 
these high excitation lines and at the same time not overproduce 
file 1 19 and 163 jjm OH lines. 

HD 100546 (Sturm et al. 2010), a Herbig Be star without 
an envelope, was chosen for early observations to be a source 
analogous to DK Cha, but further evolved. A comparison is use- 
ful because similar species are seen in its spectrum. Emission 
of H2O and OH from the upper layers of a flared, but quies- 
cent disk, similar to that around HD 100546 JSfiirm e t al. 201C|, 
could account from some of the lines from DK Cha. However, 
accounting for the difference in distance, the CO lines from DK 
Cha are more than 10 times stronger, and higher excitation lines 
from H2O are seen toward DK Cha. A possible explanation for 
the increased emission of DK Cha is an accretion shock onto the 
disk, where high densities and temperatures are present. Another 
possibility is a shock in the walls of the outflow. This will also 
be further examined in Green et al. in prep. 

4. Conclusions 

Continuum models that fit the data are consistent with the inter- 
pretation of DK Cha as transitional between a Stage I and Stage 
II object, with nearly equal masses in envelope and disk and a 
nearly face-on geometry. The full PACS SED mode observations 
reveal an extremely rich spectrum of over 50 ionic, atomic, and 
molecular lines. With further improvements to the spectral re- 
sponse function, we will be able to identify additional gas, ice 
and dust features in the spectrum. The full range scan/SED mode 
can probe many aspects of a protostellar object, from the ambient 
envelope to gas excited by ultraviolet radiation and by shocks. 
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